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ABSTRACT 

A nonlinear kinetic theory of cosmic ray (CR) acceleration in supernova rem- 
nants is employed to investigate the properties of the remnant SN 1987A. It is 
shown that a large downstream magnetic field Bd ~ 10 mG is required to fit 
the existing observational data. Such a strong field together with the strong 
shock modification due to CR backreaction provides the steep and concave ra- 
dioemission spectrum and considerable synchrotron cooling of high energy elec- 
trons which diminish their X-ray synchrotron flux below the observed Chandra 
flux which has to be considered as an upper limit for nonthermal X-ray emis- 
sion. The expected 7-ray energy flux at TeV-energies at the current epoch is 
2 x 1CT 13 erg/(cm 2 s). 

Subject headings: (ISM:) cosmic rays - acceleration of particles - shock waves - 
(stars:) supernovae: individual (SN 1987A) - radiation mechanisms: non-thermal 
- X-rays: individual (SN 1987A) 

1. Introduction 

Supernova (SN) 1987A occurred in the Large Magellanic Cloud. It has been extensively 
studied in all wavelengths from radio to 7-ray. The initial short outburst of radio emission 
(Turtle et al. 1987) is attributed to the synchrotron emission of electrons accelerated by the 
SN shock propagated in the free wind of presupernova star, which was the blue supergiant 
(BSG) (Chevalier & Fransson 1987). After about 3 years radio emission was detected again 
(Staveley-Smith et al. 1992; Gaensler et al. 1997) as well as the monotonically increased X- 
ray emission (Gorenstein et al. 1994; Hasinger et al. 1996). This second increase of emission 
is attributed to the entrance of the outer SN shock into the thermalized BSG wind and then 
in the H II region occupied by much more dense matter, consists of the swept up wind of 
red supergiant (RSG) progenitor star (Chevalier & Dwarkadas 1995). 
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To describe the observed properties of nonthermal emission of SN 1987A we use here a 
nonlinear kinetic theory. It couples the particle acceleration process with the hydrodynamics 
of the thermal gas (Berezhko et al. 1996; Berezhko & Volk 2000). Therefore in spherically 
symmetric approach it is able to quite definitely predict the evolution of gas density, pressure, 
mass velocity, together with the energy spectrum and the spatial distribution of cosmic ray 
(CR) nuclei and electrons at any given evolutionary epoch t, and the properties of the 
nonthermal radiation produced in SNRs due to these accelerated CRs. The application of 
this theory to individual SNRs (see Berezhko 2005, for a review) has provided the explanation 
of the observed SNR properties under the assumption of strong interior magnetic field. It was 
also shown that this strong field leads to the concentration of the highest-energy electrons in 
a very thin shell just behind the shock. Recent observations with the Chandra and XMM- 
Newton X-ray telescopes in space make it possible to resolve such filamentary structures, 
which are the result of strong synchrotron losses of the emitting multi-TeV electrons in strong 
magnetic fields downstream of the outer SN shock (Berezhko 2005). 

Application of the nonlinear theory to SN 1987A (Berezhko & Ksenofontov 2000) ex- 
plained the properties of the observed radioemission. It predicted also the continuous increase 
of radioemission with roughly constant rate during the period from 1999 to 2006 yr, which 
was confirmed by further measurements (Manchester et al. 2002). 

We analyze here radio data for more extended period of observation (Manchester & 
Gaensler 2001; Manchester et al. 2002) together with the Chandra measurements of X- 
ray emission (Park et al. 2004) which provides a strong constraint on the amount of the 
nonthermal synchrotron emission at the same energies 0.5-10 keV. It is demonstrated that 
this constrain gives the evidence of strongly enhanced (amplified) magnetic field inside the 
SNR, that in turn provides a better fit of radio data than in our previous study. 

2. Physical parameters of SNR 1987A 

We use canonical values of stellar ejecta mass M e j = 1OM , distance d — 50 kpc, 
hydrodynamic explosion energy E sn = 1.5 x 10 51 erg (e.g. McCray 1993). During an initial 
period the shell material has a broad distribution in velocity v. The fastest part of this 
ejecta distribution can be described by a power law dM ej /dv oc v 2 ~ k . We use a value 
k = 8.6 appropriate for SN 1987A (McCray 1993). The interaction of the ejecta with the 
circumstellar medium (CSM) creates a strong shock there which heats the thermal gas and 
accelerates particles diffusively to a nonthermal CR component of comparable energy density. 

Strongly nonuniform CSM is a result of interaction of progenitors winds: dense slow RSG 
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wind and subsequent fast and diluted BSG wind. During initial period after the explosion 
the SN shock propagated in the free BSG wind and reached after about a day ~ 1000 the 
termination shock in BSG wind, situated at the radial distance = 3.1 x 10 17 cm (Berezhko 
& Ksenofontov 2000). After that the SN shock propagates in the thermalized BSG wind of 
density p R - Considerable rapid decrease of the SN shock speed occurred during the days 
1500 — 2500 (Gaensler et al. 1997) showed that SN shock entered the H II region of density 
Pr 3> Pb- Therefore we model the CSM density distribution at distances r > R T in the form 

Pb + Pr Pb - Pr , r - R c n , 
Po = ^ 2 l c ' ^' 

where Rc = 6.1 x 10 17 cm is the contact discontinuity between two winds, lc = 0.05Rc is 
the scale of the smooth transition between them. 

The analysis of the rapidly increasing X-ray emission provides the evidence that H II 
region in turn is not uniform: at day ~ 6200 the SN shock begins to interact with the dense 
inner ring. Its radial behavior can be represented in the form 

Pm + PR2 Pry - Pm . ,r - R R . 
PR = tanh - , (2) 

where R R = 6.8 x 10 17 cm, l R = 0.08-Rr, p R2 ~ 20p R1 . Such a gas radial profile is close 
to what was extracted from the above analysis (Park et al. 2004). We adopt here for the 
gas number density N g = p/m p the following values: Ng = 0.29 cur 3 , N^ 1 = 280 cm~ 3 
and N^ 2 = 4000 cm -3 . Here m p is the mass of proton. The value N^ 1 = 280 cm which 
is by a factor of 1.5 lower than was used in our previous study (Berezhko & Ksenofontov 
2000), provides a good compromise between the SN shock dynamics seen in radio and in 
X-ray emissions (see below). Note, that the actual structure of the H II region, especially the 
dense inner ring, is essentially nonspherically symmetric. Nevertheless since our spherically 
symmetric theory reproduces the SN shock evolution R s (t) consistent with the existing data, 
one can expect that it gives a good estimate for the production of CRs and their subsequent 
nonthermal emission. 

A rather high downstream magnetic field strength Bd ~ 1 mG is needed to reproduce the 
observed steep radio spectrum (Berezhko & Ksenofontov 2000). We believe that the required 
strength of the magnetic field have to be attributed to nonlinear field amplification at the 
SN shock by CR acceleration itself. According to plasma physical considerations (Lucek & 
Bell 2000; Bell 2004), the existing CSM magnetic field can indeed be significantly amplified 
at a strong shock by CR streaming instabilities. In fact, for all the thoroughly studied young 
SNRs, the ratio of magnetic field energy density 5g/87r in the upstream region of the shock 
precursor to the CR pressure P c is about the same (Volk et al. 2005). Here B = B d /a is 
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the far upstream field presumably amplified by CRs of highest energy, a is the total shock 
compression ratio. Within an error of about 50 percent we have Bq/(8ttP c ) k5x 10~ 3 . CR 
pressure in young SNRs has a typical value P c 0.5p Vg, therefore we adopt here upstream 
magnetic field 

B = v/2tt x 10-VoK 2 - (3) 

Since the process of magnetic field amplification is not included in our theory we simply 
postulate the existence of far upstream field B Q given by Eq.(3). The spectrum of CRs 
produced by strong modified shock is very hard so that CRs with highest energies have a 
largest contribution in their energy content. These the most energetic CRs produce field 
amplification on their spacial scale that is the precursor size. Therefore CRs with lower 
energies already 'see' the amplified field B . We will argue that such a high field, that is by 
a factor of ten larger than used in our previous study (Berezhko & Ksenofontov 2000), is 
indeed required to fit the X-ray data. It also gives a better fit of radio data. 

We start our consideration from the SNR evolutionary epoch t = 1000 d, when the 
outer SN shock has a radius Ri = Rt and speed Vi = 28000 km/s. These values of R s and 
V s according to our calculations (Berezhko & Ksenofontov 2000) correspond to the end of 
SN shock propagation in the free BSG wind region r < Rt- We neglect the contribution of 
CRs accelerated in the region r < Rt, because due to a high gas density the number of CRs 
produced in the region r > Rt very soon becomes dominant. 



3. Results and Discussion 

Calculated shock radius R s and speed V s shown in Fig. la as a function of time are in 
satisfactory agreement with the values obtained on the basis of radio and X-ray measure- 
ments. Note that radio data compared with X-ray data gives larger shock size at any given 
time t and calculated radius R s (t) goes between these two sets 

To fit the spectral shape of the observed radio emission we assume a proton injection 
rate rj = 3 x 10~ 3 , which is a fraction of gas particles involved into the acceleration at the 
from of SN shock. This leads to a significant nonlinear modification of the shock: as it is 
seen in Fig. lb total shock compression ratio a 5.3 is essentially larger and a subshock 
compression ratio a s 2.8 is lower than classical value 4. 

Since the SN 1987A is very young CRs despite of very efficient acceleration accumulated 
only 7% of the explosion energy E sn (see Fig.lc). Note that the injection is expected to be 
strongly suppressed at the quasiperpendicular part of the shock, therefore one should renor- 
malize the results for the nucleonic spectrum, calculated within the spherically symmetric 
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model. The lack of symmetry in the actual SNR can be approximately taken into account by 
a renormalization factor f re pa 0.2, which diminishes the nucleonic CR production efficiency, 
calculated in the spherical model, and all effects associated with it (see Volk et al. 2003, for 
details). With this renormalization the CRs inside SN 1987A already contain 

E c = 0.07 f re E sn w 2 x 10 49 erg. (4) 

Since during the last 15 yrs the shock speed and the gas density are almost constant 
CR energy content grows roughly linearly with time. It gives the natural explanation of the 
linear increase of the radio emission detected during this evolutionary period. 

Strongly modified SN shock generates CR spectrum N oc p~ 7 , which is very steep 
at momenta p < m p c, with index 7 = (a s + 2)/(a s — 1) ~ 2.7. CR electrons with such 
a spectrum produces synchrotron radioemission spectrum S v oc v~ a with spectral index 
a = (7— 1)/2 pa 0.9, that very well corresponds to the experiment, as it is seen in Fig.2, where 
we present synchrotron energy spectra vS v , calculated for five subsequent epoch together 
with the experimental data. Note that CR spectrum has a concave shape: it becomes flatter 
at higher momenta p. As a consequence synchrotron spectrum S u {y) is also concave as it is 
clearly seen in Fig.2 at v < 10 12 Hz. Radio data reveal this property in good consistency 
with theoretical prediction. 

Strong downstream magnetic field pa 15 mG, that corresponds to the upstream field 
B pa 3 mG (see Fig.l), provides synchrotron cooling of electrons with momenta p > 10m p c 
that in turn makes synchrotron spectrum at high frequencies v > 10 12 Hz very steep (see 
Fig.2). Concave shape of electrons continuously produced at the shock front together with 
their synchrotron cooling lead to a formation of two peaks in synchrotron energy spectrum 
vS v . The first one at v pa 10 12 Hz corresponds to CR electron momentum p pa 10m p c 
above which synchrotron energy looses are relevant, whereas the second peak at v pa 10 18 Hz 
corresponds to the maximum momentum p pa 10 A m p c of accelerated electrons. Under this 
condition calculated synchrotron flux at frequency v ~ 10 17 Hz, which corresponds to the 
photon energy e 7 = 0.5 keV, is below the measured flux at the epochs t > 3000 d. Since the 
contribution of the nonthermal radiation in the observed X-ray emission of SN 1987A is not 
very well known (e.g. Michael et al. 2002), the observed X-ray flux has to be considered as 
the upper limit for the expected nonthermal emission. At early epoch t < 2500 d however 
the calculated flux exceeds the measured one (see the curve, corresponding t = 1970 d in 
Fig.2). This can be considered as indication, that the actual magnetic field B is few times 
larger then given by the Eq.(3). 

To illustrate the situation expected at considerably lower magnetic field we present in 
Fig. 3 synchrotron energy spectra calculated at the same set of parameters as before except 
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magnetic field, which was taken B = 200 /iG independent of time. Since the interior 
magnetic field ~ 2 mG (cr w 10) is essentially lower in this case, synchrotron losses of 
high energy CR electrons are considerably smaller compared with the previous case. Due to 
this fact synchrotron spectra considerably exceeds at any given epoch the measured Chandra 
flux. Therefore we can conclude, that the actual interior magnetic field strength is not lower 
than 5 mG. 

Projected radial profile of the nonthermal X-ray emission at high interior magnetic field 
Bd = 10 mG is characterized by an extremely sharp peak of thickness L ~ 10~ 3 R S just behind 
the SN shock. Since the thermal X-rays have much wider radial distribution it provides 
the possibility for experimental determination of the nonthermal X-ray contribution in the 
observed emission. Experimental determination of the X-ray spectrum at v — 10 18 — 10 19 Hz 
(e 7 = 5 — 50 keV) is the another possibility to discriminate the nonthermal emission. Such 
kind of measurements can be done with space instrument Suzaku. 

Completely different possibility to have the synchrotron flux, which fit the radio data 
and goes below Chandra data, is low field scenario, when the magnetic field is as low as 
Bq < 2 /iG. In such a case the cutoff frequency of synchrotron spectrum u max is lower than 
10 17 Hz and calculated fluxes go down exponentially at v > u rnax below Chandra points. 
However the value Bq < 2 /iG is unrealistically small for such a dense medium, which we 
have in SN 1987A. In addition the effective energy of electrons 

e e ~ 5V0V1 GHz)/(B d /10 /iG) GeV, (5) 

which radiate at v ~ 1 GHz in the interior field Bd = 10 /iG is about 5 GeV. The electron 
spectrum at GeV-energies is characterized by a power low index 7 « 2 that is considerably 
smaller than what is required for the observed radio spectrum. Note that as it is seen in 
Fig. 3 even the field B^ R3 2 mG is too small to have a good fit of radio data as in Fig.2. 
Therefore low field scenario should be rejected. 

Calculated 7-ray integral flux shown in Fig. 4 at all energies is dominated by the 7r°-decay 
component. Since the SN shock is strongly modified 7-ray spectrum at energies e 7 > 0.1 TeV 
is very hard: F 1 oc e~ 09 . At the current epoch the expected 7-ray energy flux at TeV-energies 
is about e 7 F 7 « 2 x 10~ 13 erg/(cm 2 s) and during the next four years it expect to grow by 
a factor of two. This flux is by a factor of five lower compared with our previous prediction 
(Berezhko & Ksenofontov 2000), because we use here renormalization factor f re = 0.2, which 
reduces the number of accelerated CRs calculated within spherically symmetric approach. 
The existence of strongly asymmetric CSM structure, which is dense inner ring, makes our 
prediction of 7-ray flux uncertain. According to the rough estimate this uncertainty is not 
very large, about a factor of two, due to the stronger SN shock deceleration in denser medium. 
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At the moment there are only upper limits of TeV emission obtained by CANGAROO 
(Enomoto et al. 2003) and HESS (Rowell et al. 2004) instruments (see Fig.4). 

4. Summary 

The kinetic nonlinear model for CR acceleration in SNRs has been in detail applied to 
SN 1987A, in order to compare its results with observational properties. We find that quite 
a reasonable consistency with most of the observational data can be achieved. 

The evidence of the efficient CR production leading to a strong shock modification 
comes from radio data. Significant shock modification leads to the steep and concave CR 
spectrum which very well fits the observed radioemission spectrum, under the condition of 
extremely high downstream magnetic field strength B d ~ 1 mG. To be consistent also with 
the Chandra measurements of the X-ray flux, which has to be considered as an upper limit 
for the nonthermal X-ray emission, even larger interior field « 10 mG is needed. Such a 
high field provides a strong synchrotron losses of CR electrons emitting nonthermal X-rays, 
that makes the high frequency part of the synchrotron spectrum much more steeper and 
provides a consistency with the experiment. 

The expected 7r°-decay 7-ray energy flux at the current epoch is about e 7 F 7 k2x 10~ 13 
erg/(cm 2 s) at energies e 7 = 0.1 — 10 TeV. Therefore the detection of 7-ray emission at these 
energies would imply clear evidence for a hadronic origin and for a strong magnetic field 
amplification inside SN 1987A. 
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Fig. 1. — (a) Shock radius R s , shock speed V s , gas density N g and upstream magnetic field 
B at the current shock position; (b) total shock (er) and subshock (cr s ) compression ratios; 
(c) kinetic energy of ejecta E e j and accelerated CR energy content E c as functions of time. 
The dotted vertical line marks the current epoch. The observed radius of the SN shock, as 
determined by radio (Manchester et al. 2002) and X-ray measurements (Park et al. 2004), 
are shown by circles and stars respectively. The scaling values are Ri = Rt = 3.1 x 10 17 cm 
and Vi = 28000 km/s 
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Fig. 2. — Synchrotron energy spectrum of SN 1987A, calculated for the five evolutionary 
epochs. The ATCA radio (Manchester et al. 2002; Manchester & Gaensler 2001) and Chandra 
X-ray (Park et al. 2004) data for three epochs are shown as well. Higher measured fluxes 
correspond to later epoch. 
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Fig. 3. — Same as in Fig.2, but the calculations are for the case of constant upstream 
magnetic field B = 200 /zG. 
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Fig. 4. — Integral 7-ray energy flux from SN 1987A, calculated for the same five epoch as in 
Fig.2. CANGAROO (Enomoto et al. 2003) and HESS (Rowell et al. 2004) upper limits are 
shown as well. 



